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New  General  and  Complementary  Energy  Theorems,  Finite 
Strain,  Rate-Sensitive  Inelasticity,  and  Finite  Elements: 

Some  Computational  Studies. 

S.N.  Atluri*  and  H.  Murakawa** 

Center  for  the  Advancement  of  Computational  Mechanics 
School  of  Civil  Engineering 

Georgia  Institute  of  Technology,  Atlanta,  Ga.  30332. 

Abstract: 

[General  variational  theorems  for  the  rate  problems  of  rate-de¬ 
pendent  finite  strain  inelasticity,  in  terms  of  the  appropri¬ 
ate  rates  of  the  first  and  second  Piola-Kirchhof f  stress  ten¬ 
sors,  the  symmetrized  Biot-Lure 1  stress  tensor,  and  their  con¬ 
jugate  measures  of  strain-rate,  are  discussed.  Certain  new 
rate-complementary-energy  principles,  involving  the  rate  of 
spin  and  the  rate  of  the  symmetrized  Biot-Lure'  stress  tensor 
as  variables,  are  stated  for  finite  strain  analysis  of  rate- 
sensitive  materials,  such  as  those  exhibiting  elasto-visco- 
plastic  and  creep  behavior.  Uniqueness  and  stability  criteria 
for  those  inelastic  solids,  using  the  finite  element  counter¬ 
parts  of  the  new  complementary  energy  rate  principles,  are 
discussed. Computational  studies,  using  the  complementary  ener¬ 
gy  methods,  discussed  herein  include:  ( i )  bifurcation  neck¬ 
ing  and  post-buckling  analyses  of  initially  perfect  elasto- 
plastic  bars,  and  (it)  post-buckling  and  large-deformation 
analyses  of  thin  elastic  plates  under  inplane  compression  and 
transverse  bending  loads. 

Introduction: 

The  topic  of  rate  (incremental),  multi-field,  variational 
principles,  in  general,  and  the  rate  complementary  energy 
principles,  in  particular,  and  the  corresponding  finite  ele¬ 
ment  methods,  for  finite  strain  analysis  of  compressible  non¬ 
linear-elastic  solids  were  discussed  in  detail  by  Atluri  and 
Murakawa  [1],  Also  discussed  in  (1]  were  the  contributions  of 
Koiter,  Zubov,  and  Fraeijs  de  Veubeke,  dealing  with  the  sub- 
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ject  of  complementary  energy  principles,  governing  the  total 
deformations  of  semi-linear  and/or  nonlinear  compressible  iso¬ 
tropic  elastic  materials.  It  was  shown  in  [1],  that  the  con¬ 
cept  of  treating  the  angular  momentum  balance  condition  as  an 
a  posteriori  constraint  through  a  complementary  energy  princi¬ 
ple  involving  the  symmetrized  Biot-stress  (or  what  is  also  re¬ 
ferred  to  as  the  symmetrized  Lure’ -stress  or  the  Jaumann- 
stress)  as  well  as  the  orthogonal  tensor  of  rigid  rotation,  as 
variables,  as  first  introduced  by  F.  de  Veubeke,  has  certain 
fundamentally  novel  features  that  makes  it  attractive  for 
practical  application. 

The  ideas  of  discretizing  the  angular  momentum  balance  condi¬ 
tions  through  a  complementary  energy  principle  has  been  ex¬ 
tended  by  the  authors  in  (i)  the  incremental  (rate)  analysis 
of  finite  strains  in  compressible  as  well  as  incompressible 
nonlinear  elastic  materials  [2-4],  (ii)  the  rate  problems  of 
classical  (rate-independent)  finite  strain,  elastq-plasticity 
[5-9],  and  (Hi)  nonlinear  stability  and  post-bifurcation 
analysis  of  semilinear  isotropic  elastic  beams  [10]. 

In  the  present  paper  the  authors'  earlier  work,  [2-10],  is  ex¬ 
tended  to  the  cases  of  finite  strain  analyses  of  materials 
with  rate-sensitive  behavior  such  as  elasto-viscoplasticity 
and  creep,  and  post-buckling  and  large-deformation  behavior  of 
structural  members  such  as  plates  and  shells,  undergoing  large 
rotations  and  large  stretches. 

The  summary  of  the  topics  presented  in  the  following  is:  (i) 
a  discussion  of  general  (multi-field)  variational  principles, 
with  emphasis  on  complementary  energy,  in  terms  of  alternate 
stress-rates  and  conjugate  measures  of  strain-rate,  for  rate- 
sensitive  inelastic  materials,  ( ii )  rate  complementary  energy 
potentials,  for  the  chosen  stress-rates,  for  rate-dependent  as 
well  as  rate-independent  materials,  (Hi)  criteria  for  unique¬ 
ness  and  stability  of  solutions,  (iv)  numerical  study  of  neck¬ 
ing  of  an  initially  perfect  elasto-plastic  bar,  and  (y)  numer¬ 
ical  study  of  post-buckling  of  an  axially  compressed  plate  of 
semilinear  isotropic  elastic  material,  undergoing'  large  rota¬ 
tions,  as  well  that  of  a  thin  plate  undergoing  large  displace¬ 
ments  due  to  transverse  loading. 


PRELIMINARIES: 

We  use  a  fixed  rectangular  cartesian  coordinate  system,  and 
employ  the  notation:  (_)  denotes  a  second-order  tensor;  (  ) 
denotes  a  fourth-order  tensor;  (-)  implies  a  vector;  a  =  A.b 
implies  a.  =  A..b.;  A.b  implies  a  product  such  that  (A.B).. 

j-  lj  D  ~  ~  m  ~  ID 

=  A .  .  B  . .  ;  A: B  =  trace  (A  .B)  =  A . . B . . ;  and  u. t  =  u . t . . A  parti- 

cle  in  the  undeformed  body  has  a  position  vector  x  =  x  e  (a=l..3) 
where  eQ  are  unit  cartesian  bases.  The  gradient  operator  V° 
in  the  undeformed  configuration  CQ  is  V°  =  (e^/Bx^)  .  The  po¬ 
sition  vector  of  the  particle  in  the  deformed  configuration, 

N 

say  CM,  is  y  =  y.e..  The  gradient  operator  in  CM  is  v  = 
(e^3/3y.).  The  deformation  gradient  tensor  is  F  =  (v  y)  , 
such  that  F.  =  y .  h  (3y./3x  ).  The  nonsingular  F  has  the 
polar-decomposition,  F  =  a. (I+h)  where  the  rotation  a  is  or¬ 
thogonal  and  the  stretch  h  is  symmetric  and  +ve  definite.  The 

T  T  T 

Green-Lagrange  strain  is  g  =  1/2 (F  .  F-I)  =  1/2 (e+e  +e.e) 

o  T  ~  **  **  ** 

where  e  =  (V  u)  ,  u  =  y-x . 

For  the  present  purposes,  we  introduce  the  stress  measures  ( i ) 
the  "true"  Cauchy  stress  t;  ( ii )  a  weighted  tensor,  the  Kirch- 
hoff  stress  tensor  a  =  Jt  where  J  is  determinant  of  matrix 
(y^  ];  {Hi)  the  first  Piola-Kirchhof f  stress  tensor  t;  {iv) 

the  second  Piola-Kirchhof f  stress  tensor  s,  and  (y)  the  sym¬ 
metrized  Biot  stress  tensor  (or  what  is  also  often  referred  to 
as  the  symmetrized  Lure',  or  the  Jaumann  stress  tensor)  r.  As 
discussed  in  [1,5],  and  elsewhere,  the  above  stress  measures 
are  related  as: 

1  1  T  1 

t  =  jF.t  =  jF.s.F1  =  ja  (1) 

t  =  s.FT  =  J(F_1.t);  s  =  J(F“1.T.F’"T)  (2) 

r  =  i(t.a+aT.tT)  =  j[s.  (I+h)  +  (I+h)  . s]  (3) 

The  tensors  t,  a,  s,  and  r  are  symmetric,  while  t  is  unsymme- 

-T  _  l  t 

trie.  In  the  above,  F  =  (F  )  and  the  superscript  T  de¬ 
notes  a  transpose. 

RATE  FORMULATIONS: 

Now,  we  consider  the  (incremental)  rate  analysis  of  finite 
strain  problems  of  an  inelastic  solid  with  a  rate-sensitive 
constitutive  law.  In  doing  so,  one  can  choose  an  arbitrary 
reference  frame.  In  practice,  however,  two  choices,  one  the 
so-called  total-Lagrangean  (TL)  and  the  other,  the  so-called 


updated-Lagrangean  (UL)  reference  frames  are  appealing.  Even- 
though  the  choice  of  a  reference  frame  does  not,  per  se,  af¬ 
fect  the  theoretical  or  computational  approaches,  we  discuss 
the  details  of  a  UL  formulation,  since  the  rate  consitutive 
relations  of  an  inelastic  solid  depend,  naturally,  on  the  cur¬ 
rent  state  of  true  stress. 


In  the  UL  formulation,  the  solution  variables  in  the  generic 

state  CN+1  are  referred  to  the  configuration  of  the  body  in 

the  immediately  preceding  state,  CN,  which  is  known.  In  the 

UL  formulation,  one  is  essentially  concerned  with  an  initial 

stress  problem:  the  initial  "true"  stress  in  CM  is  the  Cauchy 

stress  x  ,  while  the  initial  displacements  in  C  as  referred 

N  ^ 

to  CM  are,  obviously,  zero.  Let  y.  be  the  current  spatial  co- 

W  1  N 

ordinates  of  a  particle  in  CM.  Let  V  be  the  gradient  operator 

M  »  *’ 

in  CN  (ie. ,  V  =e^3/3y*/)  and  let  u  be  rate  of  deformation  (ve¬ 
locities)  from  CN»  We  define  the  rate  of  displacement  gradi¬ 
ent  e  =  (v*5u)T  and  write  £  =  e  +  u  where  e[e.  .  =  ^■(3il./3yI?  + 

3<i./3y. )  ]  is  the  symmetric  UL  strain-rate  and  a>[£>..  =  =-(3Ci./3y. 

3  ~  13  ^  1  3 

-3Uj/3yi)  is  the  skew-symmetric  spin-rate.  Let  x,  and 

a (=JN[x+ (e : I) xN] ) (where  JN=pQ/pN,  PQ  and  pn  being  the  mass- 
densities  in  CQ  and  CN#  respectively) ,  be  the  substantial  de¬ 
rivatives  of  the  Cauchy  and  Kirchhoff  stresses  respectively. 

As  is  well-known,  these  stress-rates  are  not  objective.  Let  t, 

s,  and  r  represent  the  appropriate  stress  rates  referred  to  C  ; 

~  ~  .  N+l  N  N+l  .  N 

ie.,  for  instance,  §At  =  s^  -x  where  sN  is  the  second 

Piola-Kirchhof f  stress  in  CN+^  referred  to  (and  measured  per 

unit  area  in)  CN«  It  is  shown  in  [5]  that: 


...  N  N  .Tw-N  1  ..  .  N.  /tN 

s  =  (g-e.g  -o  .e  ) / J  ;  t  =  (g-e.g  )/J 

N  .  ,.T  N.  .  1 /  N  ...  N. 

r  =  ^-(t+t  +x  .oj+u  .x  )  =  s  +  2 '(l  »E+e.x  ) 


(4a, b) 
(5a, b) 


Unless  large  elastic  deformations  of  a  dilatational  nature 

have  preceded  the  inelastic  straining,  one  may,  without  sig- 

N 

nificant  error,  asstime  that  J  =1.0. 

The  equations  of  linear  momentum  balance  (LMB) ,  angular  momen¬ 
tum  balance  (AMB) ,  compatibility,  and  traction  and  displace¬ 
ment  boundary  conditions  (TBC  and  DBC)  in  the  UL  rate  formu¬ 
lation  can  be  written  as: 

LMB:  VN.[§  +  xN.(VNu)]  +  pNB  =  0  (or)  VN.t  +  pNB  =  0  (6a, b) 


(7a, b) 


AMB:  s  =  sT;  (or) 

or,  equivalently, 


( VNu) T. 


N 


+  t  =  t  +  T 


N 


(VNu) 


.  N  .  N  •  r T 

to.T  +  e .  t  +  t  =  t 

compatibility: 

,  N  . 
+  t  .  e 

N  .T 

+  T  .0) 

(7c) 

e  =  e+  w=(Vu)  ; 

(or)  | 

,  nN . \  .  , _N. .  T 

=  (7  u)  +  (V  u) 

(8a, b) 

TBC: 

n*.  [s  +  tN. (VNji) ]  = 

n*.[tj 

'  • 

=  t  at  S 
-  o  N 

(9) 

DBC: 

• 

u  =  u  at  S  „ 

-  uN 

(10) 

Let  us  suppose,  for  the  moment,  that  the  consitutive  law  for  a 
rate-dependent  material  can  be  expressed  (as  shown  later  in 
this  paper)  in  terms  of  certain  rate-potentials,  as: 


s  =  3W/3e;  t  =  3U/3feT;  t  =  3Q/3e  (lla-c) 

We  consider  the  Legendre  (contact)  transformations  of  the  type: 

s : e  -  ft(e)  =  S*(§);  tT:6-£j(e)  =  E*(t) ; 

i:i  -  0(e)  =  R*(£)  (12a-c) 

such  that 

3S*/3s  =  e;  3E*/3t  =  6T;  3R*/3f  =  e 


As  discussed  in  [1-6],  the  AMB  conditions  are  embedded  in  the 
structure  of  W  and  0.  As  shown  in  [5],  the  complementary  en¬ 
ergy  principles,  and  the  Hellinger-Reissner  type  principles, 
involving  the  stress-rates  S,  t,  and  r  are  as  below.  In  each 
case  the  functional  whose  stationary  condition  is  the  prin¬ 
ciple  in  question  is  given.  The  respective  functionals  are, 
denoted  by  tt  with  the  subscript  C  and  HR  denoting  complemen¬ 
tary  and  Hellinger-Reissner  type  functionals,  respectively. 


irc(u,§)  =  J  {-£>*(!)  -  ^-xN:  [  (VNu)  .  (VNu) T]  )dv  +  j  t.  uds 

*%  SUN  (13) 

tthr(u,|)  =  J  {  -S*  ( s)  +  pNB.u  +  |tN:  [  (VNu)  .  (VNu)T] 

VN 

+  Js:[(vNu)  +  (VNu)T]}dv  -  /  t.uds  -  J  t. (u-u)ds  (14) 

SoN  SuN 


it  (t)  =  /  -£*(t)dv  +  f  t.u 

~  "  ~  “'s 


ds 


(15) 


V 


N 


UN 


xRR(t,u)  =  j  {  -E* ( t ) dv  -  pNB.u  +  tT: [ (7Nu)T] }dv 
VN 


+  contd. 


-  f  t.uds  -  f  t.(&-u)ds  (16) 

“  ^un" 

11^(4, u)  =  J  {-R*(r)  +  jtN:  (“T*^)  ~  tT:u}dv  +  J"  t.uds  (17) 

VN  SuN 

11  hr  ^ '  4 '  — )  =  J"  {“&*(£)  +  j=rJN:  (w^.w)  -  p^B.O  +  t**:  [  ^ul^-t*  :<L}dv 
VN 

-  /  t.Ms  -  f  t.(Ci-u)ds  (18) 


In  the  above  VN,  SaN,  and  SuN  are  the  volume,  prescribed-trac¬ 
tion  boundary,  and  prescribed  displacement  boundary,  respec- 
tively,  of  the  solid  in  CN;  and  p  and  B,  are  respectively, 
the  mass-density  and  rate-of-body-force  in  CN-  The  above 
functionals  are  valid,  in  general,  for  non-conservative  (de- 
formation-dependent)  surface  tractions.  If  A  is  the  set  of 
conditions  that  are  satisfied  a  priori,  and  B  is  the  set  of 
those  conditions  that  are  satisfied  a  posteriori  in  the  vari¬ 
ational  principle,  for  each  of  the  above  functionals  the  sets 
A  and  B  are  as  follows:  ( i )  Eq.  (13):  Set  A (Eqs .  6a,  7a,  and 
9),  set  B (Eqs.  8b,  and  10)  ( ii )  Eq.  (14):  Set  A  (the  exist¬ 
ence  of  S*  such  that  3S*/3|  =  e ,- and  Eq.  7a),  set  B (Eqs .  6a, 
8b,  9,  and  10).  (iii)  Eq.  (17):  Set  A (Eqs .  6b,  9,  the  de¬ 
finition  of  r  as  in  Eq.  5a  and  that  u  is  skew- symmetric) ,  set 
B (Eqs .  7c,  8a,  and  10);  (iv)  Eq.  18:  Set  A(the  definition  of 
r  as  in  Eq.  (5a),  and  that  y  is  skew  symmetric),  set  B(6b,  7c, 
8a,  9 ,  and  10) . 


The  complementary,  and  Hellinger-Reissner  type  principles  as 
through  Eqs.  (17)  and  (18)  were  first  stated  by  Atluri  [5]. 

The  general  invalidity  of  the  principles  through  Eq.  (15) 
which  was  alluded  to  by  Hill  [11],  and  Eq.  (16),  were  discus¬ 
sed  in  [5].  Eventhough,  Eqs.  (13)  and  (14),  and  the  attendant 
variational  principles,  may  be  viewed  as  being  consistent,  the 
limitations  of  practical  applicability  of  these  are  discussed 
in  [5].  Especially,  in  the  application  of  Eq.  (13),  the  need 
to  select  a  symmetric  s,  such  that  Eq.  (6a)  (which  involved 
coupling  with  u)  is  satisfied,  a  priori,  is  not  an  altogether 
easy  proposition.  Several  interesting  ways  of  satisfying  Eq. 
(6a),  and  of  application  of  Eq.  (13),  were  discussed  by 
Atluri  [6,7,12]. 


However,  the  complementary  principle  of  Eq.  (17) ,  introduced 
in  [5],  has  several  attractive  features  for  practical  appli¬ 
cation:  ( i )  the  LMB,  Eq.  (6b)  can  be  easily  satisfied  by  set¬ 
ting:  t  =  VNx£+t.P  where  \j>  are  first-order  (once  differenti¬ 

able)  stress  functions  ( ii )  u>  can  be  selected  to  be  skew  sym¬ 
metric,  by  setting  di^  =  where  e^j  is  an  alternating 

tensor).  In  general,  even  in  a  TL  formulation,  the  constraint 
T 

a. a  =  I  is  easily  met  by  taking  a  to  be  a  function  of  the  3 
Euler-angles  of  rigid  rotation  [4].  In  the  case  of  plates  and 
shells,  the  concepts  of  a  finite-rotation  vector,  as  discussed 
later,  may  be  employed. 

The  application  of  the  complementary  energy  principle  as 
through  Eq.  (17) ,  and  its  TL  rate  counterpart,  is  illustrated 
later  in  this  paper. 


RATE  POTENTIALS  FOR  RATE- SENSITIVE  MATERIALS: 

As  discussed  in  [5],  and  elsewhere,  the  principle  of  object¬ 
ivity  is  met,  in  writing  the  rate  constitutive  law  of  the  ma¬ 
terial,  by  postulating  the  constitutive  relation  between  the 
objective  strain-rate  §  and  the  objective  co-rotational  (or 
also  at  times  referred  to  as  the  Zaremba,  or  the  rigid-body  or 
the  Jaumann)  rate  of  Kirchhoff  stress,  denoted  here  by  o*.  It 
is  well-known  that. 


.  N  N  .T 

o*  =  cr-w.a  -  a.&J 

Thus,  in  view  of  Eqs.  (4-5), 

....  N  N  . .  ,  N  • 
s  =  (a*-e.0  -0  .  e)  /J  ;  t 

r  =  (o*— =-(5.0  +o  .  e )  ] / J 


(19) 


,.  .  .  N  N  .  w  N 

(o*-e  .  o  -o  .  u>)/J 


(21) 


Thus,  if  V  is  the  postulated  rate  potential  for  5*  such  that 
3V/3  e  =  g*,  we  can  define: 

JNW  =  V  -  aN:(e.e);  3W/3e  =  s  (22) 

JNU  =  W  +  joN:(eT.e);  3U/3eT  =  t  (23) 

JN0  =  V  -  |oN:(e.e);  3Q/3c  =  r  (24) 

From  these,  one  can  establish  S*(s),  E*(t)  and  R*(r)  as  de¬ 
fined  in  Eqs.  (12a-c) .  Thus,  we  focus  attention  on  the  poten¬ 
tial  V. 


It  is  worth  noting  that  for  materials  with  rate-independent 
constitutive  laws,  such  as  classical  elastic-plastic  materials 


A 


the  derivative  (')  is  considered  to  be  with  respect  to  a  fic- 
tition  time.  However,  for  rate-sensitive  materials,  such  as 
elasto-viscoplastic  and  creeping  materials,  the  derivative  (‘) 
is  w.r.t.  to  natural  time. 


For  rate-independent  classical  elastic-plastic  materials.  Hill 
[13]  presented  the  postulation: 


V  -  -L, 


2  ijk£ Eij ek l 


where  ^  i-s  a  +ve  definite  symmetric  (under  ij 


(25) 

-  k l  inter¬ 
change)  tensor  of  instantaneous  elastic  modulii,  a  =  1  or  0  ac¬ 
cording  to  whether  A^  is  positive  or  negative,  g  is  a 

scalar  related  to  the  measure  of  hardening,  and  A„  is  a  tensor 
normal  to  the  interface  between  elastic  and  plastic  domain  in 
the  space.  Prandtl-Reuss  type  rate  equations  of  type  (25) 
for  classical  isotropically  hardening  materials  can  easily  be 
derived,  formally,  to  be  [5]: 


o  (e • a' ) a' 

g*  -  2 ye  +  A ( e : I)  I  -  12au  ~  ~  ~ _  (26) 

(o’: a’) [6y+2(3Fo/3WP) ] 

where,  X  and  y  are  Lame1  constants,  o’  =  g-l/3(a:I)I  is  devi- 
atoric  Kirchhoff  stress,  and  the  yield-surface  is  represented 
by  F  =  [3J2(g’)  ]**  -  Fq  =  0. 

Here  Fo=Fq(Wp)  where  Wp  =  /g:ePdt;  and  J2(g’)  =  (1/2)  (a '  :  a ' )  . 

A  rate-sensitive  constitutive  law  of  a  considerable  generality, 
as  given  by  Perzyna  [14],  can  be  easily  written  for  finite 
strains,  when  an  associative  flow-rule  is  used,  as: 

ea  =  Y<  4>  (F)  >~  (27) 

where  <  <|>>  denotes  a  specific  function,  such  that  <  4>>  =  <)>  (F) 
for  F> 0 ,  and  <j)=0  for  F<_0.  The  parameter  y  is  called  the 
fluidity  parameter  and  e  is  the  general  anelastic  strain.  It 
has  been  shown  by  Zienkiewicz  and  Coworkers  [15],  and  Argyris 
and  Coworkers  [16]  that  classical,  rate-independent  elasto- 
plastic  solutions  can  be  obtained  from  the  above  theory,  when 
(i)  either  y->«°  or  ( ii )  a  stationary  solution  of  the  visco¬ 
plastic  flow  is  sought.  Various  forms  of  4>  were  reviewed 
by  Perzyna  [ 14) .  For  the  Hencky-Mises-Huber  yield  criterion, 
one  can  define  F  to  be: 


F 


[3J2(a')  ]** 


F  h  a 

o  eq 


F 

o 


(28) 


where  aeq  is  the  equivalent  Kirchhoff  stress.  A  simple  choice 
for  $ (F)  can  be: 

■MF)  =  Fn  (29) 

With  Eq.  (29) ,  the  viscoplastic  strain  rate  as  in  Eq.  (27)  can 
be  seen  to  correspond  to  the  well-known  Norton’s  power  law  for 
steady-state  creep  when  Fq->-0.  We  now  derive  rate  potentials 
AV (=VAt) , AW ,  etc.  for  the  viscoplastic  constitutive  laws  given 
by  Eqs.  (27-29) . 

.  N 

Let  at  times  t_,  and  tXT+At,  the  Kirchhoff  stresses  be  a  +aAt  = 

N  N  N  ~  ~ 

a  +Aa,  respectively,  where  a  is  the  substantial  derivative. 

The  inelastic  strain-rates  corresponding  to  Eqs.  (27-29)  at 
times  t^  and  t^+At,  are  given,  respectively,  by: 

%  5  <Sa>N  =  YFNOF/3g)  =  Y(3/2aeg)  (aeq-Fo)ng’  (30) 

N 

In  the  above  a'  =  (a  ) '  (ie.,  superscript  N  dropped  for  con- 
2 

venience)  and  a  =  (3/2) a': a*.  Likewise  at  a  time  At  later, 

eq  '  - 

§N+1  =  y(3/2)  (aeq+Aaeq)~1(aeq+A0eq‘Fo)n(?^A5,)  (31) 

By  straight  forward  algebra,  it  can  be  shown  that, 

e*  =  c^+Y(3/2a  ){-(Aa  /a  )  (a  -F  )na’+(a  -F  )nAa' 
~N+1  ~N  '  '  eq'  eq  eq  eq  o  -  v  eq  o 


+  n(a  -F  )n-1a’Aa  } 
eq  o  -  eq 


However,  it  can  be  easily  shown  that 

Aa  =  (3/2) [ (a ' : Aa ' ) / a  ]  (33) 

Using  (33)  in  (32)  we  obtain: 

Sn+!  -  Sn  +  Y:-A?  <34> 

wherein  the  definition  of  V  is  apparent.  Since  A a  can  in¬ 
clude  the  effects  of  pure  spin  between  t^  and  tN+^,  one  can 
replace  Eq.  (34)  by: 

Sn+i  "  +  Y:A?*  (35) 

where  Aa*  =  a*At  is  the  corotational  increment  of  Kirchhoff 
stress.  Now,  the  corotational  rate  g*  can  be  written  as: 


g*  =  Le: (E-e ) 

Where  L  is  the  tensor  of  instantaneous  elastic  moduli. 
=  e 

t  +At 

Thus,  Aa*  =  L  :\c  -  L  :  f  c  dt 

=  e  e  i-J  ~ 

N 

one  may  use  the  approximation. 


when  Eqs.  (35)  and  (38)  are  used,  Eq.  (37)  becomes. 


Aa  *  =  Le:Ae  -  A tLe : ( e^+BVAa*)  (39) 

From  which,  upon  rearranging  terms, 

Ac*=M:(Ae-Ae)  (40) 

wherein  the  definition  of  M  is  apparent,  and  Ae  is  known  and 

a  3  ~ a 

is  given  by:  Ae^  =  AteN-  From  Eq.  (40)  one  can  immediately 
write 


AV  =  ^ijk^'ij^kl 


^ijk£AeaijAek2. 


(41) 


From  Eq.  (41) ,  the  potentials  AW,  AU,  and  AQ  can  easily  be  ob¬ 
tained  through  Eqs.  (22-24)  . 


It  is  noted  that  Wang  [17]  attempted  to  derive  a  relation  bet¬ 
ween  As,  Ae,  Ae  ,  and  hence  AW  directly.  However,  this  deriv- 

*'  **  **  cl 

ation  appears  to  be  in  error,  since,  among  other  reasons,  the 
transformation  between  the  deviatoric  part  of  s  and  the  devi- 
atoric  part  of  t  was  assumed  to  be  the  same  as  that  between  s 
and  t  themselves. 


All  the  above  developments  for  the  UL  rate  formulation  can  be 
converted  to  a  TL  rate  formulation  by  noting  the  relations  [5]: 


E’  =  (FN)T.e.FN;e'  =  e.FN;  s'  =  JN(FN)  1.S.(FN)“T 

t'  =  JN(FN)_1.t;  r'  =  (1/2) [tN.a’+a’T.tNT+t' .aN+aNT.t,T] (42) 


where  E'  is  TL  rate  of  Gree-strain,  and  e',  s',  t',  and  r'  are 

TL  rates  of  e,  s,  t,  -and  r  respectively.  Now  a’  is  subject  to 

NT 

the  constraint  that  a  .oc  is  skew-symmetric. 


UNIQUENESS  &  STABILITY  CRITERION: 

In  the  present  paper,  the  application  of  the  rate  complementary 
energy  principle  as  embodied  in  Eq.  (17),  or  its  TL  counter¬ 
part,  will  be  used  in  some  computational  studies.  In  the  di¬ 
rect  application  of  Eq.  (17),  the  assumed  stress  field  t  must 
not  only  satisfy  the  LMB  condition  within  each  element,  but 
also  the  traction  reciprocity  condition  at  the  interelement 
boundary,  viz.,  ( n* . t) ++ (n* . t)  =  0  at  p  (where  +  and  -, 
respectively,  indicate  the  two  sides  of  ,  the  interface 

between  mth  and  (m+l)th  elements  in  C^) .  In  the  present  work 
this  interelement  condition  is  introduced  as  a  posteriori  con¬ 
straint,  through  a  Lagrange  multiplier  u  at  in  the 


functional  of  Eq.  (17),  thus  leading  to  a  'hybrid'  finite  ele¬ 
ment  method.  The  thus  modified  functional  is: 


irHQ(w,t,&  )=E{  f  [-R*(r)  +  (l/2)xN:  (wT.-l)-tT.a)]dv 
Hb  ~  -p  m  J  ~  ~ 

-  ■  mN 

+  f  (n*.t).uds  +  f  (n*.t) .u  ds  (43) 

umN  mN 

At  the  point  of  bifurcation,  or  instability,  from  the  concept 
of  adjacent  compatible  states,  the  following  criterion  can  be 
shown  to  hold: 


with 

^HS  =  0  and  5itHS 

the  constraints: 

„N  •  _  .T 

V  ,t  =  0;  w  =  -w  ; 

0 

n  *  .  t  =  0  at  S  ; 

amN 

(44) 

and 

u  =  0  at  S  „ 
umN 

(45) 

In  the  case  of  linear  pre-buckling  states,  the  above  criterion 
reduces  to  an  eigen-value  problem,  with  the  eigen-value  depen- 
ding  on  t  . 

EXAMPLE  PROBLEMS: 

Necking  of  an  Initially  Perfect,  Plane-Strain,  Elastic- 
Plastic  Bar. 

Cowper  and  Onat  [18]  examined  the  above  bifurcation  necking 
problem,  for  a  bar  of  rigid-plastic  work-hardening  material, 
under  uniform  tension  applied  at  the  ends  of  the  bar.  Mises' 
yield,  and  isotropic  hardening  criteria  were  used  [18] .  In 
[18]  only  the  eigen-value  problem  for  the  applied  tension  at 
which  necking  would  initiate  in  the  bar  was  treated,  but  the 
phenomenon  of  post-bifurcation  necking  was  not  treated  in  [18]. 
Lo(L)  and  Bo(B)  are  the  initial  (current)  length  and  width  of 
the  bar  respectively.  If  y^  and  y^  are  the  current  cartesian 
coordinates  of  a  material  particle,  the  ooundary  conditions 
are:  (i)  at  y^=0:  u  =-V;  t  =0;  [ii)  at  y^=L:  u  =+V;  t  =0; 

N  1  .!•*•  i  ^  11/ 

( Hi )  at  y2=+ (B/2)  :  t^=t2=0. 

In  the  present  analysis,  the  complementary  energy  formulation 
based  on  the  TL  rate  equivalent  of  Eq.  (17)  was  used.  The 
problem  parameters  used  are:  (B  /L  )  =  (1/3);  t  (yield 

4  o  o  y 

stress)  =  4x10  psi;  the  true  stress  versus  logarithmic  strain 

(In (1/1  ))  curve  was  assumed  to  be  bilinear,  with  the  two 
°  7  4 

slopes,  E=10  psi,  and  h  =  5x10  psi.  The  notation 


0 


=  (L-L  )L  ;  is  used.  The  eigen-value  solution  for  n  at  bi- 
o  o 

furcation-necking  for  the  perfect  bar,  as  obtained  in  [18] 
for  the  present  linear-hardening  but  rigid-plastic  material, 
is  nC  =  0.48. 
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The  variation  of  total  applied  load  with  n  is  shown  in  Fig. 

c 

1,  from  which  it  is  seen  that  necking  starts  at  n  =  0.482 
(which  is  in  excellent  agreement  with  the  value  of  n  =0.48 
of  [18]).  It  is  also  seen  that  necking  starts  after  the 
maximum  load  is  attained.  The  convergence  of  nC  with  the 
finite  element  mesh  has  been  reported  elsewhere  [8,9],  with 
the  mesh  as  shown  in  Fig.  2,  which  is  used  to  obtain  the 
remainder  of  the  reported  results,  being  the  finest  mesh 
reported  in  [8].  The  variation  of  the  width  reduction  ra¬ 
tio,  (B-Bq)/bo  is  shown  in  Fig.  2,  from  which  it  is  seen  that 
at  n  =  n°  the  width  reduction  becomes  much  more  pronounced 
as  compared  to  the  rigid-plastic  fundamental  solution. 

The  variation  of  6/L^  (with  6  being  defined  as  the  difference 
of  widths  at  loading  edge  and  the  necking. sections,  respec¬ 
tively)  with  n  is  shown  in  Fig.  3.  The  slope  of  this  curve 
at  the  beginning  of  necking,  viz,  n  =  nC  =  0.48,  was  obtained 
in  an  asymptotic  analysis  in  [18].  The  present  result  for 
this  initial  slope  is  in  excellent  agreement  with  [18].  How¬ 
ever  as  necking  develops,  the  slope  (Dn/35)  decreases  from 
the  initial  value  at  n  =  nC,  which  appears  to  be  in  contra- 


diction  with  the  result  of  McMeeking  and  Rice  [19].  The 
necked  profile  of  the  bar  for  n  >_  nC,  are  shown  in  Fig.  4. 


Figure  3  Figure  4 

The  necked  profiles  of  the  bar,  and  the  progressive  develop¬ 
ment  of  unloaded  regions  (shaded)  are  shown  in  Fig.  5  at 
various  values  of  n.  Note  that  unloading  begins  at  the  cen¬ 
ter  of  the  loaded  face  of  the  bar  at  n  =  nC  =  0.482.  Final¬ 


ly  the  distribution  of  Cauchy  stresses,  T^finthe  direc¬ 


tion  of  loading)  x22  (in  the  width  loading) , 


and  t33  (in 


the  thickness  direction  of  this  plane-strain  specimen) ,  at 
the  neck  (y^  =  L/2)  are  shown  in  Fig.  6.  These  results  are 


in  excellent  qualitative  agreement  with  those  of  Needleman 


[20]  who  also  analyses  the  necking  and  post-necking  problem 


of  an  initially  perfect  cylindrical  bar.  It  is  noted  that 


Figure  5 


Figure  6 


the  problem  of  necking  of  an  elastic-plastic  bar,  with 
initial  imperfections,  was  also  analysed  by  Osias  [21], 
McMeeking  and  Rice  [19],  and  Nemat-Nasser  and  Taya  [21]. 

It  is  observed  that  the  finite  element  meshes  used  in  [19- 
21]  are  very  much  finer  than  the  one  used  presently.  Even- 
though  a  precise  mathematical  statement  as  to  this  appears 
impossible,  the  above  comparison  of  the  mesh-sizes  appears 
to  indicate  the  relative  advantages  of  the  present  comple¬ 
mentary  energy  method. 

Finite  Deformation  and  Post-Buckling  Analyses  of  a  Thin 
Plate; 

We  consider  large  deformations  (large  rotations  and  large 
stretches)  of  a  thin  plate  made  of  a  semi-linear  isotropic 
material  (ie.,  a  material  exhibiting  a  linear  relation  be¬ 
tween  the  stretch  tensor  h  and  its  conjugate  stress-measure, 
the  symmetrized  Biot-Lure'  or  Jaumann  tensor  r) .  We  invoke 
the  well-known  Kirchhof f-Love  type  plausible  deformation 
hypotheses  that  a  normal  to  the  midplane  of  the  initially 
flat  undeformed  plate  remains  normal  to  the  deformed  mid¬ 
plane  and  that  there  is  no  thickness  stretch.  In  order  to 
derive  a  consistent  complementary  energy  principle  for 
this  constrained  deformation  problem,  we  start  with  the 
general  (Hu-Washizu  type)  variational  principle  involving 
h,  a,  u,  and  t  as  variables. 

By  introducing  the  appropriate  approximations  to  all  these 
variables,  we  derive  a  general  variational  principle  for 
the  plate  problem;  from  this  we  proceed  to  construct  a  con¬ 
sistent  complementary  energy  principle  for  the  plate  prob¬ 
lem.  It  is  shown  in  [1]  that  a  general  functional,  for 
the  three-dimensional  finite  elasticity,  whose  stationary 
conditions  lead  to  all  the  appropriate  field  equations,  is 
given  by; 


7T  (u,h,a,t)  =  /  { W(h)  -  p  cj.u  +  tT;  [  (I+Vu) 

G  -  -  ~  -  Jv  ~  o - ~  ~  -- 


-  a. (I+h) ] }dv  -  /  t.(u-u)ds  -  f  t.u 


uds 


uo 


ao 


(46) 
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where  W(h)  is  the  strain  energy  density  per  unit  of  initial 
volume  V  ,  as  a  function  of  pure  stretch  h;pQ  is  mass  density/ 
unit  initial  volume,  and  g  are  body  forces/unit  mass. 

Let  x^,  x2  be  cartesian  coordinates  in  the  mid  plane,  and  x^ 
normal  to  the  midplane,  of  the  plate.  The  position  vector  of 
an  arbitrary  material  point  in  the  undeformed  plate  is  x=x.e. 
(i=i,...3).  Under  the  present  deformation  hypotheses,  the 
position  vector  of  the  same  particle  after  deformation  is, 

Y  =  (x^+ujOe^  +  (x2+u*)e2  +  u|e3  +  x3N  where  N  is  a  unit  nor¬ 
mal  to  the  deformed  midplane.  Further,  the  displacement  u* 
(i=l,2,3)  are  functions  of  and  x2  only.  Thus  the  displace¬ 
ment  of  an  arbitrary  particle  in  the  plate,  is  y-*  =  u|e^  + 
u|§2  +  u^§3  +  x3(N-§3).  The  base  vector  at  an  arbitrary  point 


in  the  deformed  plate  are  given  by: 

ay/ax  =  (6  ;  +u*  )e.+N,  x=G  (<*=1,2;  i=l,2,3)  (47) 

3y/3x3  =  N 

The  deformation  gradient  is 

F  =  (vy)T  =  G  e  +  Ne,  (u=l,2)  (48a) 

-  —  —a— a  — — j 

Further,  for  the  present  kinematic  hypotheses,  we  assume  the 
stretch  tensor  to  be: 

^  ha&®a®6;  haB  =  *WXi}  Ii=1'***3'*  6=1,2  3  (48b) 

From  Eqs .  (48a, b)  it  is  seen: 

N  =  F.e3  =  [a. (I+h)  ].e3  =  2.e3  (49) 

Thus  the  displacement  vector  can  be  written  as: 

u(X;L)  =  u*(xa)ei+(a-I)  .e3x3  [i=l,2,3;  a=l,2]  (50) 

Further,  we  assume  that  h  Q(x.)  can  be  approximated  as: 

up  l 

hag(Xi)  =  h*aB(Xa}  +X3XaB(X6)  [a'6'Y=1'2]  (51) 

ie.  h  =  h*  +  x3 x 


For  the  semilinear  isotropic  material  we  assume  the  constitu¬ 
tive  lav/: 


r  =  2„h  +  x(h: I) I 


(52) 


Since,  for  isotropy,  h,  g,  and  r  are  coaxial,  Eq.  (3)  becomes, 
r  =  t.ct  (53) 

The  tensors  t  and  a  are  assumed  to  be: 

?  =  tijSiSji  [i. 3=1. 2, 3;]  (54) 

where,  ^  (Xj.Xj.Xj)  i  -  •ij(x1.x2)  (55) 

Finally,  the  external  forces  distributed  on  the  plate  are  as¬ 
sumed  to  be  specified  per  unit  area  on  the  mid  plane  of  the 
plate  to  be  =  g^ (x^ ) (i=l, . . 3 ,  a=l,2).  When  the  assumptions 
in  Eqs.  (50-55)  are  substituted  in  Eq.  (4-6) ,  we  find  through 
atraight-forward  algebra,  that 


*GtUi-  Ki'  * 


^6  ' 


a  . 


T  .  , 
a  i 


Mei> 


where 


and 


=  f  { W* (h* , x )  -  g-u*  +  ^T:[eaea  +  u*aea 

s 

o 

+  (a.e3)e3  -  a.(I+h)]  +  AT:[(a.e3)  aea  -  a.^Hds 

-  J*  [T.  (u*-u* )  +  M.<  (a-a)  .e3>]dc  -  C  [T.u*+ii.  (a-I)  ,e3ldc 

ter\ 

^  =  T  .< 

01' 

T  .  =  /  t  .dx,;  M  .  -  I  t  .  x_dx_. 
oi  /  ai  3  cti  /  ,  oi  3  3 

w*  =  f 

J  V 


'a  (56) 

(a=l , 2 ;  i=l,...3),  and. 


v3 
Wdx. 


(57) 


(56) ,  sQ  is  the  area  of  the  undeformed  midplane,  and 


In  Eq 

C  an 
u 

daries  of  s.  It  is  seen  that  only  t 


Cu  and  Cq  are  the  displacement  and  traction  prescribed  boun- 


ai 


enter  into  the  above 


energy  expression  due  to  the  presently  invoked  deformation  as¬ 
sumptions.  The  constitutive  equations  and  LMB  conditions  ob¬ 
tainable  from  Eq.  (57)  are: 

3W*/3h*  =  4. a  R  and  3W*/3x  =  A. a  =  N  (58) 

and  T  .  +  g.  =  0  (a=l,2;  i=l,...3)  (59) 

cu,a  J  l 

When  Eqs.  (58,  59)  and  the  appropriate  traction  boundary 
conditions  on  T  .  are  satisfied  a  priori,  one  can  eliminate 

ci-i 

from  Eq.  (56),  (i)  h*  and  x  through  the  usual  contact  trans¬ 
formations  and  by  establishing  a  complementary  energy  density 
R*  such  that  3R*/DR  =  h*  and  3R*/3N  =  x,  ( ii )  u*  through 
satisfying  (53)  a  priori.  When  this  is  done,  we  obtain  a 


complementary  energy  functional: 

^(a,  R,  N)  =  J- R*(R,N)  +  ^T:  [  eaey+ (a  .e^)  §3—3;  ] 

+  M^: [ (g . e,)  e  ]]ds-  f  (T.u*+H.<  (a-a  ).e,>)dc 
— o  , a-a  j  —  —  —  ~  —  i 


U 


-fj 


M.  (a -1 )  .  ^dc . 


(60) 


In  Eq.  (60) ,  a  is  required  to  be  orthogonal  and  further  a  is,  as 
assumed  in  Eq.  (55),  a  function  only  of  x^  and  x2»  Also  the 
variables  R  and  N  in  Eq.  (60)  are  assumed  to  be  defined  in 
terms  of  &,  and  a  as  in  Eq.  (58) .  To  assume  an  orthogonal 
a(x^,x2),  the  concept  of  a  finite  rotation  vector  [23]  is  use¬ 
ful.  Let  oj  be  the  finite  angle  of  rotation  around  an  arbitra¬ 
rily  oriented  unit  vector  e  in  the  midplane  of  the  plate.  The 
finite  rotation  vector  is  defined  to  be: 

0  =  (sinw)e  (61) 

The  action  of  finite  rotation  on  a  vector  V  can  be  written 
[23]  as  the  transformation  of  V  to  V*  as. 


V*  =  V  +  Q  x  V  +  [2x(fixV]/2cos  (oj/2)  =  g.V 
where  a=I+2xJ+  [ (£xl ) • (ftxl ) ] /2cos2 (w/2) 

It  can  be  shown  that  a  of  Eq.  (63)  is  orthogonal,  ie. , 
T 

a-a  =  I. 


(62) 

(63) 


The  vector  e  in  Eq.  (61)  can  be  written  as: 

e  =  e^cos0  +  e2  sine  (64) 

Thus,  the  rotation  tensor  a  of  Eq.  (63)  is  a  function  of  two 

parameters:  0(x^,x2)  and  u(x^,x2).  The  explicit  expression 

for  a  can  be  shown  to  be: 

2 

a=<  1-  ( 1-cosw)  sm  6>e^e^+<  ( 1-cosw)  sinecosO ^^e,, 

+  (sinwsinO) e^e2  +  ( l-cosw) sin0cosOe2e^ 

2 

+  <  1- ( l-cosw)  cos  0>e2e2  -  sinwrosO^^  -  sinwsinoc^e^ 

+  sinucosoej^  +  008106363. 

In  a  Von-Karman  type  plate  theory  w  is  assumed  to  be  moderate- 

2 

ly  large,  such  that  cosu  -  l-(w  /2)  and  sm.i  =  u;  while  the 
angle  0  can  be  assumed  to  be  arbitrary. 

Further  details  of  the  analysis  of  large  rotations  and 
stretches  of  thin  plates  using  the  complementary  energy  method 
sketched  above,  which  are  omitted  here  for  space  reasons,  will 
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